function S1map =sigmoidS1(S1map,sat)
nMap = numel(S1map);
for iMap =1:nMap 
    S1map{iMap}=sat.*(2./(1+exp(-2/sat.*S1map{iMap}))-1);
end